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Abstract 

Coupled DGLAP equations involving singlet quark and gluon distribu¬ 
tions are explored by a Taylor expansion at small x as two first order partial 
differential equations in two variables : Bjorken x and t {t = In^). The 
system of equations are then solved by the Lagrange’s method and Method 
of Characteristics. We obtain the proton structure function F 2 {x,t) by com¬ 
bining the corresponding non-singlet and singlet structure functions by both 
the methods. Analytical solutions for F 2 {x, t) thus obtained are compared 
with the recent data published by HI and ZEUS as well as with NNPDF3.0 
parametrization and their compatibility is checked. Comparative analysis fa¬ 
vors the analytical solution by Lagrange’s method. 
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1 Introduction 

In our earlier work [1], we have made an extensive comparative study on the appli¬ 
cability of the two analytical methods: Lagrange’s method and method of charac¬ 
teristics in context of the unpolarised non-singlet structure function . We have 
discussed different kinematic regions in which both the methods showed validity by 
comparing with both data and exact results for non-singlet sector. These analytical 
methods can be extended simply to embrace the singlet sector, so that it can be used 
to hnd the gluon and quark distribution functions. Solutions of DGLAP [2H5] evo¬ 
lution equations give quark and gluon structure functions which produce ultimately 
proton, neutron and deuteron structure functions. The standard program to study 
the X dependence of quark and gluon PDFs is to compare the numerical solutions of 
the DGLAP equations with the data and so to £t the parameters of the x profiles of 
the PDFs at some initial scale Qg and the asymptotic scale parameter A. However, 
for analyzing exclusively the small-x region, there exists alternative simpler analysis, 
yielding analytical solutions of the DGLAP equations [6HI0]. Some approximated 
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analytical solutions of DGLAP evolution equations suitable at small-a:, have been 
reported in recent years with considerable phenomenological success. 

Following the procedure of ref. pQ, in this present communication we have solved 
the integrodifferential equation for the quark and gluon distribution functions in the 
leading order (LO), using the analytical methods and construct our analytical solu¬ 
tions for proton structure function F,[{x, t) as the sum of a flavor nonsinglet and 
a flavor singlet F^ distribution. In the experimental front, HI [13] and ZEUS [H] 
collaboration has published very recent data on F:[{x,t) covering wide kinematic 
region of Bjorken x, and low to medium four-momentum transfer squared, Q^. In 
HI [T3| collaboration data, the earlier measurements of F^^x^t) are superseded by 
the recently published data in the kinematic region 6.5 x 10“^ ^ x < 0.32 x 10“^ 
and l.bGeV'^ < < SOOGeU^, whereas in ZEUS [H] F^ix^t) is measured in the 

region 0.00025 < x < 0.00493 with 9GeU^ < < llOGeU^. In this paper we have 

discussed the relative compatibility of both the analytical methods in context of 
F 2 {x,t) at leading order (LO) with the recent data from HI and ZEUS experiment 
as well as with the latest NNPDF3.0 [T5| parametrization based on HERA Run-H 
data. A newer data-set allows the extension of the kinematic range towards lower 
values of x for comparative analysis, in which our analytical models are approxi¬ 
mated at. 

In section 2 we describe the formalism, section 3 is devoted to testing our prediction’s 
comparison with the data, while in section 4, we give our conclusion. 


2 Formalism 


2.1 Singlet coupled DGLAP equations in Taylor approxi¬ 
mated form 

The coupled DGLAP equations for quark singlet (S(x,Q^)) and gluon {G{x,Q‘^)) 
densities are ms], 

d ( S(a:,Q2) \ ^ g, (Q2) 

91n(5^ \ G{x,Q‘^) J 27r 


P P 

J- qq J- qg 

P P 

^ gq ^99 


E(x,Q2) 

G (x, Q2) 


( 1 ) 


where as (Q^) is the strong coupling constant, PjjS are the Altarelli-Parisi splitting 
functions and the symbol 0 stands for the usual Mellin convolution in the hrst 
variable dehned as, 

a{x) ® f{x) = [ —a{y)f . (2) 

Jx y \yJ 

Introducing the variable t = In ^ and using the explicit forms of the splitting 
functions Pij{i,j = q,g) in LO, the evolution equation for singlet distribution can 
be written as |16j . 


dFi{x,t) _ ^ 
dt t 


{3 -I- 41n(l — x)}F 2 {x^t) -I- 2 


dz 


(1 - z) 

2F^{x, f)} -f |n/ /J dz {z^ + (1 - zf) G (f, t) 


{l + z^)F^ 


= 0 


(3) 
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Similarly the DGLAP evolution equation for non-singlet distribution can be written 
as, 


dF^^{x, t) 
dt 


^ [{3 + 41n(l 


+2 




= 0 . 


(4) 


Here Af = ^,/3o = being the number of flavors considered and a gif) = 

t) and F^^{x, t) are the singlet and non-singlet structure functions of the 
proton. We write Eq.(|3]) and Eq.(jl]) as, 


~ ^ [{3+ 41ti(l + ^(1^.*)] = 0 (6) 

= A [{3 + 4 ln(l - X)} F™(x. i) + ir(x, i)] (6) 

where 



To carry out the integrations in Eqs.([7}l9|), we introduce the variable u dehned as 
u = 1 — z and expand the argument — as a series. 


x 

z 


x 

1 — u 


= X 


= X+X ^ 


u 


k=0 


k=l 


( 10 ) 


Since x<z<l,so0<n<l — a;; hence the series is convergent for |m| < 1 and we 
can use Taylor expansion of t), c^nd F^^{^,t) ^ approximated form 

and as x is small in our region of discussion, the terms containing and higher 
powers of x can be neglected as those terms are still smaller and therefore, we can 
rewrite. 





kdF.^’^^i 


u 


X, t) 


k=l 


dx 


( 11 ) 


G(^.t)^G(x.t)+x^rF-FM. ( 12 ) 

z ox 

k=l 

where the terms containing x^ and higher powers of x are neglected at small x. 
Using the above Eq. flTT]) and Eq. flT^ we carry out the integrations in 2 ; in Eq.([7}l9]). 
Neglecting terms 0{x‘^) which is justihed at small x , we get. 


jf (x, t) 'zz (2x 


3)F2^(x,f) -h 


^x -|- 2x In 



dFi{xA) 

dx 


(13) 
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(l - G{x,t) - 
2 ln(—) + (1 — 

X 
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5x — 3x In — 

X 


dG{x, t) 
dx 


— (x — 1) (x + 3) F^^{x, t), 


(14) 

(15) 


The exact relation between the gluon distribution G{x,t) = xg{x,t) and quark 
distribution F^^x^t) = {qi{x,t) + gj(x,t)} is not derivable in QCD even in 

LO. However, simple forms of such relation are available in literature to facilitate 
the analytical solution of coupled DGLAP equations. In ref Cl. it was assumed 
that dependence of both the distributions are identical. In ref. [18] , on the other 
hand, the following simple relation was assumed. 


G{x,Q‘^) = k.Fi{x,Q‘‘) 


(16) 


where parameter k has to be determined from experiments. 

Again as the input singlet and gluon parametrization, taken from global analysis of 
parton distribution functions, which incorporate different high precision data, are 
also functions of x at hxed Q^, hence relation between singlet structure function 
and gluon parton densities may be expressed as a function of x [T^ . 

However a more rigorous analysis was done by Lopez and Yndurain m and they 
investigated the behavior of the singlet F^i^x^Q'^) and gluon G{x,Q‘^) as x —>■ 0. 
They observed that 

Fi{x,Q\^o = Bs{Q^)x-^^ (17) 

G{x,Q%^o = Bg{Q‘^)x-^^ (18) 

where Bs and Bq are dependent as Xq = Xs and Xs is strictly positive. Thus, 


G{x,Q^) 


fm 


It suggests a more general form [21] following. 


G{x,Q^) ^ K(Q^)Fi(x,Q^-) 


(19) 

( 20 ) 


than Eq. llTHil . 

Now using Eq. (lT3il and Eq. llTTjl we can express Eq.(l3]) in a more precise form as. 


dFi{x,t) _ ^ 

dt t 


3 + 41n(l — x)F 2 *(x, t) + (2x — 3)^2* (x, t) + x + 2xln — 

X 


1\ 9^2^(x,f) 


t 


Uf{l-lx)G(x.t)-^i5x-3xln 


dx 

1\ dG{x,t) 
x / dx 


Similarly using Eq. flTS]) we write for the non-singlet F^^ as. 


dF^^{x,t) AfX dF^^{x,t) 


dt 


dx 


21n(—) + (1 — x^) 

X 


^ [3 -I- 41n(l — x) 

NS/ 


+ {x-l){x + 3)]F^^{x,t) = 0 ( 22 ) 


0 ( 21 ) 
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Eq.(j22]) is a partial differential equation for the non-singlet structure function t) 

with respect to the variables x and t. The solutions of the Eq. fl2^ by both Lagrange’s 
and Method of Characteristics has been reported in our earlier work [T]. So we don’t 
discuss the solution of the non-singlet structure function F^^{x,t) any further in 
this work. We continue our discussion to obtain the analytical solutions for Eq. (l2T]) . 
Using above relation given by Eq. fl20|) we express Eq.([2T]) as, 


dFi{x,t) 

dt 


^ {3 -I- 41n(l — x) -I- (2x — 3)} F^^x, t) + Uf ^ K(Q‘^)F2 {x, t) 

dFi{x,t) 


t 


2x\F--”-L 

X 2 


3xln- K{Q^ 

X 


dx 


= 0(23) 


which is a partial differential equation for the singlet structure function F 2 {x,t) 
with respect to the variables x and t. We solve this PDE Eq. fl2^ with the two 
formalisms described here, the Lagrange’s method and Method of Characteristics. 
In order to do that we express Eq. fl2^ as. 




dt 


dx 


where 


cui = + 2x\n — — — ( 5x — 3xln — ) KtQ"^)} 

3po X 2 \ X J 

W 2 ^ + 41n(l - 2 ^) + {2x - 3) + n/ (l - ^j-) K{<^)) 


(24) 

(25) 

(26) 


2.2 Solution by the Lagrange’s Auxiliary Method 

To solve the equation Eq. fl24p by the Lagrange’s Auxiliary method [22], we write 
the equation in the form. 



Q(xX^^%^'‘Kp(xF^^f'^^^R(x.t,Fi{x.t)) 

(27) 

where 

Q{x, t) =t 

(28) 


F(x, t) = —ui 

(29) 

and 

i?(x, t, F^{x, t)) = R'{x)F^{x, t) 

(30) 

with 

R'{x) = UJ2 

(31) 


The general solution of the Eq. (l27l) is obtained by solving the following auxiliary 
system of ordinary differential equations. 


dx dt dF^^x^t) 

P{x) Q{t) R{x,t, F2{x,t)) 
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If u{x, t, F 2 ) = Cl and v{x, t, F 2 ) = C 2 are the two independent solntions of Eq.(l27]), 
then in general, the solntion of Eg . 0271) is 

F{u, v) = 0 (33) 

where F is an arbitrary fnnction of u and v. 

In this approach we try to hnd a specihc solntion that satishes some physical con¬ 
ditions on the strnctnre fnnction. Snch a solntion can be extracted from the combi¬ 
nation of u and v linear in F 2 , the simplest possibility being, 

u + av = f3 (34) 


where a and are two qnantities to be determined from the bonndary conditions 
on F 2 . Solving Eq. fl32l) . we obtain. 


and 


u{x, t, F 2 ) = tX^{x) 

(35) 

v{x,t,Fi) = Fi{x,t)Y^{x) 

(36) 

and Y^{x) are dehned as 



(37) 


(38) 


The explicit analytical form of X^{x) in the leading (i) approximation (at very 
small X region log(-) a:log(-) 3> x) comes ont to be. 


X^(x) = exp[ 


6(3o 


A{A + 3nfK{Q^ 


■log [log x]] 


(39) 


Using the physically plausible boundary conditions for structure functions,i.e. 

F 2 {x,t) = ^ 2 ^( 3 ;, to)i for t = to (40) 

F 2 {l,t) = 0 for any t. (41) 


and putting the values of u and v in Eq. llMll . we obtain the solution for Eq. (1T711 as. 


F^{x,t) = F^{x,to) (J- 

, ix»(x) - xyi)] 

(42) 

1 |.YS(i) - (i).YS(l)| 

Eq.(|39j) gives us. 



X»(l) 

= 0 

(43) 

which yields. 



F^{x,t) = F^{ 

(44) 


Eq. (l44ll gives the t evolutions of singlet structure function at LO and is the solution 
for F 2 by Lagrange’s method. 
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2.3 Solution by the method of characteristics 


To solve the PDE Eq. (j2T|l by the method of characteristics [23l|2l], we express it in 
terms of a new set of coordinates (s,r), such that Eq. fl2n) becomes an ODE w.r.t. 
one of the new variables. We know that most of the important properties of the 
solution of Eq. fl2T|) depends on the principal part of the equation i.e. the left hand 
side in Eq. fl27p . This part is actually a total derivative along the solution of the 
characteristic equation, 


dx ui 

dt t 


(45) 


which gives the characteristic curves of Eq. fl24p . That is along the characteristic 
curve, the partial differential equation becomes an ordinary differential equation. 
The characteristic equation Eq. fldSp can be written as. 


dx dx ds 
dt ds dt 

with, 

dt 

ds 


(46) 

(47) 


dx 

ds 


= -Ui 


(48) 


Using Eq. fHHl) in Eq. fl2Tl) . the left hand side becomes an ordinary derivative with 
respect to s and the equation becomes an ordinary differential equation. 


where 


dFi{s,T) 

ds 


+ {s,t) F^{s,t) = 0 


(49) 


C^{s,t)=UJ2 (50) 

0 J 2 is as given above in Eq.([26]). 

As per the initial conditions i.e. x{s = 0) = r and t{s = 0) = to, the solutions of 
the characteristics equations yield. 


with. 


s = loglog(-)“i 

X 

t .A.' 

T = X exp I — I “1 




0 


CKi — 


3/3o 


4{2 + A-(02)|} 


(51) 


(52) 


Integrating Eq. (H^ along the characteristic curve, we obtain the solution for F^i^x, t) 
in (s, r) space as. 



( 53 ) 
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where, 


.^1 = 4 log ( 1 — rexp 


' ^ ^ — 

- 

'to 


+2r exp 


to 


+% ( 1- 


- 

'to^ 




(54) 

Using the solutions of the characteristic equations Eq. flST]) . which lead us to the 
{x,t) space, we can express Eq. flS^ in a more precise form as. 


F/(j,i) = Ff(i,t„)(f)"W> (55) 

to 

where, 

(56) 

and F^^t) = F^^Xjto) is the input function. Eq. fl55|l is the analytical solution for 
the singlet structure function within the present formalism. Unlike Eq. fl44|) . Eq.(|55]) 
is sensitive to gluon distribution as well as KlQ"^), which occurs in spite of the 
leading log approximations (at very small x region (log(-) S> a:log(-) 3> x)). The 
reason is that in the Lagrange’s method, F.^^ [1] and F 2 has identical evolution 
for the least approximated level. It is possible only when the gluon effect is absent. 
But this feature has already been well observed [I8l[25], the only new observation is 
that it is true even when k = K{Q‘^), i.e. dependent. 

Using our results derived in this section, we will calculate the proton structure 
function F^{x, t) from the singlet and the non-singlet structure function using the 
relation 

Fi = —F^^ + —^2^ (57) 

2 18 2 18 2 K J 

We discuss in the next section the phenomenological consequences of our results 
derived in this section. 


3 Form of K(Q‘^) 

Let us now discuss the plausible forms of K{Q‘^) as dehned in Eq. fl20|) above and 
discuss the related constraints on the parameter duos k and a. 

3.1 Choice of the form of K{Q‘^) 

The important characteristics of pQCD is the logQ^ dependence, as can be seen 
from the dehnition of running coupling constant, as well as any evolution of 
structure function. However alone does not appear and this basically yields the 
following plausible forms for the function iL((5^), 

1. 

= (58) 

which is a power law in Q^. 
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II. 


AW) = fc(log|l 


2\ O' 


(69) 


•which gives us logarithmic dependence in Q^. 


III. 


K{Q^) = k{ log log — 


(60) 


allowing us log log dependence. 

These forms are related to the other two possibilities through the correspondence, 


log t log log t 
jj = a —^— = v- 


t 


logt 


(61) 


While dehning the forms of it is reasonable to have a growth of ‘f’ instead 

of Q^. Possible generalized forms of K{Q‘^) are, 


K{Q^) = kJ2Ci 


i=l 


9 : 

Ql 


2\ Mi 


K{Q^) = kJ2Q{ty 


i=l 


K{Q^) = kY,Q [log 


2=1 


(62) 

(63) 

(64) 


We have chosen the function K{Q‘^) with logarithmic dependence in i.e. Eq.( 
for dehniteness and simplicity. But such a pattern comes with large number of 
parameters. However such proliferation of parameters makes the phenomenology 
uninteresting. Hence economy of parameters in terms of numbers justify that the 
most appropriate and QCD inspired functional form for the function K{Q‘^) has to 
be of the logarithmic form and we consider it to be. 


A'(02) = fc(log|))"=H' 

where k and a are two parameters to be hxed. 


(65) 


3.2 Reality constraint on the parameters k and a 

The essential condition for our analytical solution for F^, obtained by method of 
characteristics, Eq. fl55]l to be real is that the exponent n{x,t) has to be real. This 
imposes a reality condition on which quantitatively leads us to the conclusion 
that. 


0 < rexp 



< 1 


( 66 ) 
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which also follows from the dehnition of r as given by Eq. (l5Tll . The new variable r 
is dependent on both k and a as ai is parametrized by them. The inversely propor¬ 
tionate feature of cti on both k and a makes the realization that k and a cannot be 
too large. So for any value of x and the choice of our parameter duos k and a is 
bounded by the above reality condition and cannot be treated as free parameters. 
The recent HERA data allows us to explore a wide range of evolution,i.e. 
1.5 < < SOOGeV^ for F^^x^t). Hence given the two parameters k and a, 

we £x them for this entire region of Q^. The reality condition allows an effective 
range of values for both k and a for the considered region, the best htted range 
for k and a being 0.001 < h < 1.45 and 0.001 < cr < 0.055 respectively. 


4 Results and Discussion 

In this particular work we have calculated the evolution for singlet structure 
functions using two analytical approaches. The proton structure function F^{x,t) 
has been calculated using the relation Eq. flFTIl . extending our comparative study 
of the above discussed analytical methods in terms of proton structure function 
F^{x,t) with very recent experimental data published by HI [13] and ZEUS [H] 
collaboration. We have used the LO MSTW 2008 [26] input for evolution of our 
solutions with Qo = IGeV^. We have plotted our analytical solution for f-evolution 
of F^{x,t) by Lagrange’s method with HI [I3] and ZEUS [H] experimental data in 
hgure 1 and hgure 2 respectively. We have considered the range 10““^ "G x < 10“^ 
and 2.5GeU^ < < 500GeU^ for HI data and 0.00025 < x < 0.00493 and 

9 < < llOGeU^ for ZEUS data. As in case of ZEUS data, the extracted values 

of F 2 have been given only for 27 (x, bins against few hxed values, hence 
we have considered the entire data range while plotting. Good consistency has been 
observed between our analytical solution by Lagrange’s method and the experimen¬ 
tal measurements by both HI and ZEUS within the range 0.00025 < x < 0.001 and 
3 GeV^ <Q^< Q 0 GeV\ 

In figure 3 we have shown the evolution of our analytical solution for F^{x, t) by 
method of characteristics, along with the HI experimental data. As can be seen the 
evolution of the analytical solution by method of characteristics for F^{x,t) is not 
compatible with data, further we have observed that structure function has been 
decreasing with increasing Q^. This behavior has been observed to be true for other 
regions of x too as well as for ZEUS data, hence not included in the text here. 

In hgure 4 we have plotted the analytical models with the very recent NNPDF3.0 [15] 
parametrization against for different values of x. We have conhned the compari¬ 
son within the region 0.00025 < x < 0.013 and 3GeV^ < < 85GeU^, for F^ix, t) 

produced by the NNPDF3.0 collaboration based on the HERA-H data. Here the 
vertical error bars represent uncertainties given by the standard deviations and com¬ 
puted by added in quadrature method in our work. We note that comparison with 
NNPDF3.0 parametrization also supports our observation with data and it is clear 
that the region of validity of our analytical solution by Lagrange’s method is ap- 
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proximately in the range 0.0001 < x < 0.01 and 3GeV^ < < OOGeld^. However 

our analytical solution by method of characteristics for does not follow the 

general growth of evolution. 

As per pQCD prediction structure function should rise in small x with piTlI^. 
But from the form of as has been dehned in Eq. flMj) . it is evident that for x < 0.66 
it is always positive and the exponent of our solution Eq.(|55]), i.e. n(x, t) is negative. 
Thus the behavior of Eq. dSB]) contradicts both QCD expectation [HlElISZj and data. 
Further if we extrapolate (^i to large x i.e. for x > 0.66, then it is observed that 
.^1 is negative yielding a positive exponent n{x,t), which yields actually expectation 
of small X QCD. The expected large x behavior of QCD is that structure function 
should fall [27] . So the prediction extrapolated to high x is not even consistent with 
QCD prediction. We therefore infer that method of characteristics is less favored 
than Lagrange’s method as per the approximated solutions in DGLAP approach 
within the considered formalism, a feature which has presumably been overlooked 
earlier. 

We note the sensitivity of the two parameters k and a along with the Eq. fl20|) to¬ 
wards F^ix, t), which have played a crucial role in case of method of characteristics. 
However in case of the solution by Lagrange’s method, given by Eq. fHl]) . such re¬ 
striction as k and a does not appear in the solution. Effectively it means that 
Lagrange’s method allows such growth with Q^. 


5 Conclusion 


This work is an extension of the work of ref. |T] for the two important analytical 
methods, Lagrange’s and Method of Characteristics, in obtaining the analytical solu¬ 
tions for proton structure function F.[{x,t), which consists of non-singlet 
and singlet F 2 {x,t) structure functions. For this part we pursue a general form as 
given by Eq. fl20|) . relating F 2 {x,Q‘^) and G{x,Q‘^) for comparison with theoretical 
analysis of [20] . However consequence of the relation could not be tested separately 
in the present work as it does not effect the analytical solution by the Lagrange’s 
method, obtained in the leading log - approximation and only the solution by 
method of characteristics has exclusive dependence on the relation. We summarize 
our comparative analysis by looking at the impact of sensitive comparison with the 
three different sources of predictions for F 2 {x,t), HI, ZEUS and NNPDF3.0. How¬ 
ever data analyzed in the range 0.00025 < x < 0.001 and 3GeV^ < < 60GeU^ is 

found to favor the former (Lagrange’s method) and not the later (method of char¬ 
acteristics) . 

We conclude this section with a comment. We have followed the results of the gen¬ 
eral analysis of ref. [18] and used in [27] and not incorporated the observation of 


ref. [28]; i.e. for x —)■ 0 and —)■ oo, the ratio ^ ^ 


Ini^) 


1 . 

2 


, which suggests 


that the exponent a of Eq. fl65l) might have additional x and dependence as well. 
Such effects are currently under study. 
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Figure 1: Proton structure function F^{x, t) as a function of for different fixed x 
values by Lagrange’s method with HI data. Here the dash-connected line represents 
our analytical model. 
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Figure 2: Proton structure functionF 2 ^(x, t) as function of by Lagrange’s method 
with ZEUS data. Experimental data is considered for different bins of x and Q^. 
Here the line represents our analytical model. 
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Figure 3: Proton structure function as a function of for different fixed 

X values by method of characteristics with HI data. Here the dash-connected line 
represents our analytical model. 
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Figure 4: Proton structure function F^{x,t) as a function of for different hxed 
X values by Lagrange’s method and method of characteristics with NNPDF3.0 data 
with standard deviation. Here the dashed and dotted line represents our analytical 
model by Lagrange’s method and method of characteristics respectively. 


15 




















